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Résumé. — On étudie la complexité et la stabilité des calculs dans la jaco- 
bienne des courbes de grand genre sur le corps des complexes avec une attention 
particulière aux courbes modulaires. 

Abstract (Jacobians and numerical stability). — This paper is concerned 
with the complexity and stability of arithmetic opérations in the jacobian variety 
of curves over the field of complex numbers, as the genus grows to infinity. We 
focus on modular curves. 
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1. Introduction 

Il est traditionnel de calculer dans le groupe des points de la jacobienne d'une 
courbe algébrique projective lisse et géométriquement irréductible X de genre g en 
représentant tout élément de ce groupe par un diviseur effectif de degré g, une fois 
choisi un tel diviseur O comme origine. La somme de deux diviseurs P —O eX Q — O 
est réduite par le calcul de l'espace linéaire associé au diviseur P + Q — O suivi de 
la localisation des zéros d'une fonction non nulle de cet espace. 

Comme l'application de Jacobi, 



S<^X ^ Jx 

de la puissance symétrique g-ième S^X de X dans sa jacobienne Jx, n'est pas un 
isomorphisme, la représentation n'est pas unique. 

Classification mathématique par sujets (2000). — llFll, 11F25, 11F30, 11Y16, 11Y35, 
65E05, 65Y20, 68Q15. 

Mots clefs. — jacobienne, approximation, stabilité, formes modulaires, complexité algorith- 
mique, machine de Turing, temps polynomial déterministe. 
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Si le corps de base est un corps fini F,, les opérations arithmétiques y sont exactes 
et rapides. 

On considère ici le cas oti le corps de base est le corps C des complexes. On se 
donne un modèle analytique naturel et une mesure sur X{C). On s'intéresse à la 
complexité des algorithmes utilisés pour ajouter et réduire des diviseurs. Le cadre 
est celui des machines de Turing classiques. En effet, on peut avoir en vue des 
applications arithmétiques comme le calcul de nombres de points, ou de coefficients 
de formes modulaires et les calculs en nombres complexes ne sont alors qu'une étape 
dans la recherche d'une quantité discrète. Le projet de Bas Edixhoven pour répondre 
à une question de René Schoof [6, 7, 3] se prête à cette approche. 

Bien sûr, les machines de Turing ordinaires ne manipulent pas les nombres réels ni 
complexes mais plutôt des nombres rationnels, décimaux ou binaires. Cependant, on 
peut voir un nombre réel a comme un oracle qui, pour tout entier positif fc, retourne 
une valeur binaire ou décimale approchée de a à exp(— /e) près. Si une machine de 
Turing doit résoudre un problème dont les entrées sont des nombres réels, elle reçoit 
un oracle pour chacun de ces réels. Si la machine de Turing calcule un nombre réel, 
on lui donne en entrée la précision absolue k requise et elle retourne une valeur 
approchée à exp(— fc) près du résultat. On dit que la machine est polynomiale si elle 
répond en temps polynomial en la taille des données et k. On note que la recherche 
des racines complexes d'un polynôme unitaire à coefficients complexes se fait en 
temps déterministe polynomial grâce à la méthode de quadrichotomie de Weyl par 
exemple. On veut dire par là qu'une valeur approchée à exp(— fc) près de chaque 
racine peut être calculée en temps polynomial en le degré du polynôme, la taille des 
coefficients (logarithme du maximum des modules des coefficients) et la précision 
absolue k requise. 

On veut savoir si la complexité asymptotique des opérations arithmétiques dans 
la jacobienne est polynomiale en le genre de la courbe. La première difficulté est 
de donner un sens précis à cette assertion. Plutôt que de rester dans le vague, 
on formule et on étudie ces questions dans le cas important et représentatif des 
courbes modulaires Xo{p) lorsque p est un entier premier qui tend vers l'infini. 
L'algorithmique de ces courbes est riche et largement explorée. On trouve dans 
[4, 5, 8] des algorithmes pour l'étude homologique des courbes modulaires et des 
méthodes analytiques expérimentales motivées par la vérification de conjectures 
arithmétiques et la recherche de points rationnels. 

La section 2 décrit le modèle analytique standard de ces courbes ainsi que ses 
propriétés algorithmiques. On y rappelle d'abord les résultats de Manin, Shokurov, 
Cremona et Merel concernant le calcul des périodes, et on en donne une expres- 
sion quantifiée du point de vue de la complexité algorithmique et de la stabilité 
numérique. Cette dernière est assurée en dernier ressort par des minorations du 
volume des périodes et du déterminant jacobien de l'application d'intégration de 
Jacobi. Ces minorations reposent elles mêmes sur des considérations d'intégralité 
des coefficients des formes modulaires primitives, propres et normalisées. 

On présente dans la section 3 des algorithmes pour les opérations élémentaires 
dans la jacobienne Jo{p) et pour la résolution effective du problème inverse de Ja- 
cobi. La complexité et la stabilité de ces algorithmes sont étudiées avec les outils 
de la section 2 puis estimées dans les théorèmes 1 et 2. On obtient des algorithmes 
déterministes polynomiaux en p. Le caractère déterministe de ces algorithmes s'ex- 
plique en dernier lieu par la connexité du tore analytique complexe Jo(p)(C). 
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Tous les lemmes et définitions concernant la localisation et la stabilité des zéros 
de fonctions analytiques sont présentés dans l'appendice A qui est indépendant mais 
doit être au moins parcouru avant de lire les sections 2 et 3. 

Les méthodes, les énoncés et les démonstrations que nous donnons pour les 
courbes Xoip) s'étendent sans peine au cas de Xi(p). Pour les courbes modulaires 
de niveau composé, il faut une majoration des coefficients des développements de 
Fourier en toutes les pointes ainsi qu'un algorithme pour les calculer. 

On trouvera un index à la fin de cet article. 

Convention importante : le symbole O désigne partout une constante absolue 
positive et effective, chaque fois différente. La présence de ce symbole dans une 
formule ou un énoncé signifie que cette formule ou cet énoncé sont vrais si, pour 
chaque occurence, ce symbole est remplacée par une constante positive effective 
bien choisie. 



2. Courbes modulaires Xo{p) 

Cette section rappelle, précise et complète quelques résultats métriques et al- 
gorithmiques concernant les courbes modulaires Xq{p). On supposera que p est 
premier et que le genre g de Xq{p) est au moins 2. 

Le paragraphe 2.1 introduit quelques notations et un recouvrement non injectif 
de Xo{p) par deux disques analytiques centrés en chacune des deux pointes. 

Les propriétés élémentaires des formes primitives, propres et normalisées sont 
rappelées dans le paragraphe 2.2 et celles de l'homologie dans le paragraphe 2.3. 
Le calcul des périodes est abordé dans le paragraphe 2.4. Ces trois paragraphes 
résument le travail de Manin, Shokurov, Cremona et Merel sur cette question. 

Le paragraphe 2.5 établit une minoration du volume du réseau des périodes. Une 
formule d'intégration sur les surfaces de Riemann relie ce volume au produit des 
normes de Petersson des formes primitives, propres et normalisées, ces dernières 
étant faciles à minorer parce que le développement de Fourier commence par (1 + 
0{q))dq où g = exp(2i7rT) est le paramètre de Tate associé à un r du demi-plan de 
Poincaré. 

Le paragraphe 2.6 établit des majorations simples mais nécessaires des intégrales 
de Jacobi et définit V inst,abilité d'un diviseur effectif de degré g. Le paragraphe 2.7 
construit un diviseur d'instabilité assez petite. Cela revient à trouver g points dans 
le modèle canonique de Xo{p) qui ne soient proches d'aucun hyperplan. Autrement 
dit, le jacobien en ces g points n'est pas trop petit. On prend le parti (maladroit 
en pratique mais simple en théorie) de chercher les g points dans le voisinage de la 
pointe à l'infini. Le terme principal du développement du jacobien y est le wronskien. 
On le minore grâce à l'intégralité des coefficients de son développement de Fourier. 

Le paragraphe 2.9 étudie la stabilité de l'application inverse de Jacobi. Celà se 
réduit à majorer la dificrence entre cette application et sa linéarisée. 

La connaissance d'un g-uplet de points de faible instabilité, donné au paragraphe 
2.7, permet de construire au paragraphe 2.10 des sous-ensembles finis de taille mod- 
este et bien distribues dans le tore complexe. Comme les éléments de ces ensembles 
sont images par l'application de Jacobi de diviseurs connus, il sont des auxiliaires 
précieux pour la résolution approchée du problème inverse de Jacobi. Ils permettent 
de discrétiser ce problème. 
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2.1. Un modèle analytique. — Soit p un nombre premier et X = Xç,{p) la 

courbe modulaire de niveau p associée au sous groupe de congruence F = Fq {p) de 
SL2(Z). On note H le demi-plan de Poincaré et H* = H U Q U {00} . La surface de 
Riemann compacte quotient r\W* est X{C). Son genre est g = p+'^-^^-^'^s ^^g^ 

1/2 = 1+ ^^^ etf3 = l-|-^^^.Ilyaf2 points elliptiques d'ordre 2 et on note 

V2 le diviseur somme de ces points. De même il a 1/3 points elliptiques d'ordre 3 et 
on note V3 le diviseur somme de ces points. Voir [15, Propositions 1.40 et 1.43]. Le 
genre de X est compris entre et . Le quotient r\W est un ouvert de Zariski 
de X noté Y = Yo{p). On note que la largeur de la pointe 00 est 1 et la largeur de la 
pointe est p. Pour t G H on pose q = q{T) = (/^(t) = cxp(2?'7rr) et w{t) = — -^r 
et q' = q'{T) = qoir) = (?(«;(r)) = exp(^). On note P = Poo = P(t) = P(g)'ie 
point de Y associé à r et P' = Pq = P'(t) = P'{q) = P{w{t)) = W{P) = P{q') où 
W est l'involution d'Atkin-Lehner. On a le diagramme 




On pose S = 



de sorte que St = — l/r et Tt = 



Étant donnés deux réels Poo et Pq plus petits que 1 on peut se demander si 
l'union de l'image par P^ du disque ouvert D(0,R^) et de l'image par Pq de 
P>(0,Po) recouvre X(C). 

1 )''^=[o 1 

r + 1. Soit TZ le domaine fondamental usuel de SL2(Z), délimité par le cercle de 
centre et de rayon 1 et par les droites d'abscisses —1/2 et 1/2. Alors un domaine 

fondamental pour F est constitué de l'imion de TZ et des ST^TZ pour k entier de à 
p— 1. Ces derniers sont contenus dans l'image par S de l'ensemble des r = a-\-ib avec 
b> Donc leur image par w est constituée de complexes dont la partie imaginaire 
est au moins Si on choisit Rq > exp(— ^^!^) alors l'image de D{0,Ro) par Pq 
recouvre les ST^TZ pour /s de à p — 1 . 

Comme TZ est contenu dans le demi plan des parties imaginaire au moins égales 
à \/3/2 on prend Poo > exp(— tta/S) et l'image de -D(0, Poo) par Pqo recouvre TZ. 

On pose donc Pœ = 0.005 et Pq = 1 — ^ 

On a donc recouvert X{<C) par l'image de deux disques analytiques complexes 
= D(0, Poo) et Do = £)(0, Pq). 



2.2. Différentielles. — On peut maintenant calculer des espaces de formes 
différentielles sur X. On fixe donc un entier d > 1. A toute forme modulaire 
parabolique / de poids 2d sur F on associe la différentielle oj = {2mYf{dTY de 
degré d. D'après [15, Proposition 2.16] on a 



Div(a;) = Div(/) - d(0) 



„ , d^ 2d^ 

d{^) - -V2 - yP3. 



On pose donc = (d - 1)(0) + (d - l)(oo) + [|JP2 + cherche 
une base 2?^ de l'espace Ti.'^{Ad) des formes différentielles de degré d et de diviseur 
> -Ad. 
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On prend pour Vd l'ensemble des w = {2iTiY J {q){dTY = ^{dqY où f{q) 
est une forme modulaire parabolique primitivc^^', proprc^^^ et normalisée'-'^) sur 
r = ro(p) et de poids 2d. Si / est une telle forme elle admet un développement 
/ = X^jt>i OfcÇ^ avec ai = 1 et pour tout entier fc > 1 on montre que le coefficient 
ttfe est un entier algébrique majoré en module par fc''"'"^. Il suffit de le montrer 
pour k — une puissance d'un premier £. D'après le théorème de Deligne on a 
\ai\ < 21'^-^ et d'après [1, Theorem 3] 

n(2d-l) 

donc \ain\ < Unl 2 ou u„ est la suite récurrente uq = 1, Ui = 2 et Un+2 = 

2u„+i + Un. Donc Un = (l + ^/2)"+Ml-^/2)"+l |^^| < 471 < ^^^^ |^^^| < 

„ n(2d-l) 

Le développement de u en q^o est donc le développement standard, donné par 
les valeurs propres des opérateurs de Hecke. On peut calculer les coefiicients au 
comme valeurs propres des operateurs de Hecke agissant sur les symboles de Manin- 
Shokurov suivant [4, 14, 10]. Les plongcmcnts complexes des valeurs propres peu- 
vent alors être approchés en temps polynomial en p et la précision absolue requise, 
par un algorithme de recherche de racines de polynômes comme celui de Weyl. 

Le développement de to en q^ est le tire en arrière de lu par l'application Pq : 
D — {0} Y. Comme Pq est la composée de P^o et de W, le développement de ui en 
go est le développement de W{uj) en q^. Mais w est vecteur propre de W de valeur 
propre ±1. On a donc la même majoration pour les coefficients du développement 
de oj en qo. 

Lemme 1 (Meinin, Shokurov, Cremona, Merel). — // existe un algorithme 

qui pour tous p premier, d > 1 et r > 1 calcule les plongcmcnts complexes des 
r premiers coefficients de toutes les formes modulaire primitives, propres et nor- 
malisées de niveau p et poids 2d en temps polynomial en p, d, r et la précision 
absolue requise. 

2.3. L'homologie de la courbe. — La théorie de Manin [12, 14, 4] établit 
que l'homologie relative Hi{X,ptes,Z) est engendrée par les symboles modulaires. 

Un symbole est noté indifi^érement {c : d) = ^^^ia'fl- L'ensemble des 

symboles est P = ¥\Z/pZ). On rappelle que 5* = ^ ~ ^ ^ et T = ^ ^ ^ 

Pour tout 7 G SL2(Z) on note (7) le symbole {7(0), 7(00)}. 

On note B le sous-Z- module de engendré par les {c : d) + {c : d)S = {c : 
d) + i-d : c) et (c : d) + (c : d)TS+ (c : d){TSY = {c : d) + {c + d : -c) + {d:-c-d) 
où (c : d) parcourt P. 

On note Z le sous Z- module libre et saturé de Z^ engendré par les (c : 1) pour 
c 7^ et par (00) = (0 : 1) + (1 : 0). C'est le module des symboles à bord nul. La 
base formée des (c : 1) pour c 7^ et de (00) permet d'identifier Z au réseau Z^ de 
MF muni de la forme bilinéaire canonique (de matrice identité dans cette base). 



(^^Cela signifie qu'elle ne provient pas d'une forme de niveau plus petit. Se dit en Anglais newform. 
Cette condition est vide ici puisque le niveau p est premier. 
'^■'Autrement dit, / est vecteur propre des opérateurs de Hecke. 
f^^Son développement de Fourier commence par q. 
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On note que B c Z c M^. Comme le quotient Z/B = iîi(X, Z) est sans- 
torsion, le sous module B est saturé dans Z. On identifie Hi{X,Z) à la projection 
orthogonale de Z sur le M-espace vectoriel de MP orthogonal au sous-espace vectoriel 
MB engendré par B. 

Comme B est engendré par des vcctciirs de norme < \/3 et qu'il a pour dimension 
p — 2g, son volume V est un entier positif majoré par . Mais d'après [13, 

Proposition 1.2.9, Proposition 1.3.5] Hi{X,Z) est inclus dans et son volume 

est égal à 1/V. 

D'après l'inégalité d'Hcrmitc [13, Théorème II.2.1], le réseau V'^Hi{X,Z) C 

admet une base constituée de vecteurs de norme*^^-* majorée par (|) ^ 3^^"2-^. 
L'algorithme LLL [2, Theorem 2.6.2] produit en temps polynomial en p une base 
de V'^H-i{X, Z) formée de vecteurs entiers de norme < 2^'t^3't^ < 3^ . Donc ces 
vecteurs sont des combinaisons de symboles avec des coefficients majorés par cette 
même borne. On peut faire beaucoup mieux en y regardant de plus près. Au total, 
on obtient une base B de Hi {X, Z) C ZP dont les vecteurs sont des combinaisons 
de symboles à coefficients dans ^Z dont les numérateurs sont des entiers bornés 
en valeur absolue par 3^*^ . 



2.4. Les périodes. — On observe que le diviseur Ai est nul. Donc les formes 
primitives, propres et normalisées de poids 2 donnent une base I>i de l'espace des 
formes différentielles holomorphes. 

Le réseau C des périodes est construit en intégrant chaque forme de Vi le long 
des symboles comme font Tingley dans sa thèse et Cremona dans son livre [4, 

Proposition 2.10.1]. Si <? G ^o{p) avec g = ( ^ ^ ) et c > on pose yo = et 



pc d 



xi = —dyo et X2 = ayo- Alors pour toute forme / = J2k>i poids 2 on a 

(1) / fiq) — = — exp(— 27rnyo)(exp(27rma:2) — exp(27rma;i)). 

Cette quantité est évaluée en temps polynomial en p, c et la précision absolue 
requise et elle est majorée en module par un polynôme en p et c. 

Si (c : 1) = ^ ^ ^ = {0, -1} est un symbole différent de (0 : 1) et (1:0), on 

peut supposer que 1 <c<p. Afin d'utiliser la formule d'intégration (1), on choisit 

deux entiers uetv tels que uc — vp=let on note que la matrice g = \ ^ ^ | € 

\ pv c J 

Tq{p) vérifie ^(0) = 1/c. On peut choisir Q <u,v <p. Donc /(^-i) f{l)^ est évaluée 
en temps polynomial en p et la précision absolue requise, et elle est majorée par un 
polynôme en p. 

Comme les éléments de la base B de Hi {X, Z) construite ci-dessus sont des com- 
binaisons linéaires des symboles (c : 1) et (oo) G B avec des coefficients dans ^Z 
à numérateurs majorés par 3^ en valeur absolue, on peut calculer les périodes uj 
pour 7 e ;B et u; G Pi en temps polynomial en p et la précision absolue requise et 
ces périodes sont majorées en module par exp(p^). 



Noter que dans le livre de Martinet, la norme d'un vecteur est définie comme la valeur de la 
forme quadratique en ce vecteur. Ici, la norme est la racine carrée de la forme quadratique. 
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2.5. Le tore complexe. — L'espace H = T-0'®'h} des formes harmoniques admet 
une forme bilinéaire définie par intégration. Si Wi = U\dq + Vidq et LO2 = U2dq + V2dq 
on pose < u)i,u)2 >= J-^ '^1 ^ '^2- 

L'intégration définit aussi un accouplement J : H x Hi{X,C) C qui à (0^,7) 
associe la période J^^ lu. Il en résulte un isomorphismc entre H et le dual de Hi{X, C). 
Mais l'accouplement d'intersection induit un isomorphismc entre Hi {X, C) et son 
dual. On en déduit un isomorphismc l entre H et Hi {X, C) qui à tout cj associe 
l'unique 7 = tel que J^w = -)-g pour tout g G Hi {X, C). D'après [9, Proposition 
IIL2.3.], cet isomorphismc est une isométrie : 



<LtJl,UJ2>= / LUI AlJ2 = ■ i{LU2). 

JX 

On définit l'opérateur * : H ^ H par * (udq + vdq) = {—iudq + ivdq). On définit 
un produit hermitien sur H par 

{wi,u)2)= / uiiA*Û!2 = i / {uiÛ2 + ViV2)dq A dq. 
JX Jx 

C'est le produit de Petcrsson. Les opérateurs de Hcckc sont autoadjoints pour ce 
produit hermitien. Donc deux formes distinctes dans Vi sont orthogonales. Il reste 
à évaluer (w, w) pour chaque élément uj de la base Pi. 

On note B la base de H\{X, Z) duale à gauche de B pour la forme d'intersection, 
c'est-à-dire que pour toTit 7 G S il existe un unique 7 e B tel que 7 • 7 = 1 et si 
7 9^ 7' on a 7 • 7' = 0. Ainsi = Y^^^b 7 

Donc 

{oj, oj) =< w, *û) >= i{uj) ■ l{*û)) = i{uj) ■ L{iÎL)) = i 7 • 7' / w / 



7,7'eB 

Notons P la matrice 2g x g des périodes holomorphes 



LJ. 

7 J"f' 



Soit M = {P\P) la matrice 2g x 2g des périodes harmoniques. Soit M* = {iP\ — 
iP). Soit Q = (7 • 7')^ -y'gg la matrice de la forme d'intersection dans la base B. 
Soit Q = = (7 • 7')^ la matrice de la forme d'intersection dans la base B. 

La matrice du produit scalaire de Petersson dans la base î?iU2?i de H n'est autre 
que *MQM*. Comme le déterminant de Q est 1 et comme le volume du réseau C 
des périodes est le module du déterminant de M divisé par 2^, il vient que ce volume 
est le produit des ^{oj^oj) pour w dans la base Vi. 

On peut minorer chacun des en notant que 



dq = (1 + Vafc/ 



q " ^ 

^ k>2 

avec |afc| < /c'^ de sorte que J^ujA*ûj = |1 + ^j,>2 ctfeç'^^^pîrf'Z A dq. On observe 
que le disque formé des q de module inférieur à exp(— 27r) est contenu dans un 
domaine fondamental de X donc {u},u)) > 27rr^ min|g|<r |1 + J2k>2^k'l''~^\'^ pour 
tout r < cxp(— 27r). 

Or pour |g| < r on a I Zk>i «fe+i^'l < ^^•Efe>o(fc + ^f^" < jêw ^ ^'^^^ 
r = exp(-27r). Ainsi > 27rexp(-47r)(l - 0.016)^ > 2 • lO"'"' de sorte que le 

volume du réseau C des périodes est au moins 10~^^. 
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Comme les périodes iv sont majorées par exp(p'-') on en déduit que le volume 
de tout sous-réseau du réseau des périodes est minoré par exp(— p'-'). De sorte que si 
l'on connaît les périodes avec une précision absolue polynomiale en p, on les connaît 
aussi avec une bonne précision relative. En particulier, si on connaît un point de 
C^^ par ses coordonnées complexes alors on connaît le point du tore C^^/£ puisque 
£ n'est pas trop petit ni trop aplati. 

Lemme 2 (Volume et complexité du réseau des périodes) 

Si Xo{p) est de genre g > 1 on note Vi la base de ?i^{Xo{p)) constituée des 
formes primitives, propres et normalisées et on identifie le dual de 'H^{Xo{p)) à 
C^^. On appelle réseau des périodes le réseau de C^^ formé des périodes de Xq{p). 
Ce réseau est de volum,e > 10~^^. Tous les sous-réseaux non nuls du réseau des 
périodes ont un volume > exp(— f^^^^) où ci est une constante positive effective. Le 
réseau des périodes admet une hase constituée de vecteurs de norme < exp(p'^2) où 
C2 est une constante positive effective. Une telle hase peut être calculée en temps 
polynomial en p et la précision ahsolue requise. 

2.6. L'application d'intégration de Jacobi. — On note ^,o : -D(0, Ro) ^ C^^ 
l'application d'intégration de Jacobi 



On ne précise pas l'origine o de l'intégrale. Cette application est bien définie à 
une constante additive près. On définit de même ji^o '■ D(f),R^) C''^^ en veillant 
à choisir la même origine o 

Ces intégrales se calculent par intégration terme à terme de la série associée à la 
forme différentielle. 

Le disque D(0, 1) C C est muni de la distance usuelle associée à la norme sur C 
définie par le module z t-^ \z\. 

L'espace C^^ peut être muni des normes ou L°° dont la définition est rappelée 
au paragraphe A.l de l'appendice. 

L'application d'intégration de Jacobi est Lipschitzienne dans le sens suivant. 

Soit Pi un point de coordonnées ti, qi et q[ tel que qi = exp(2i7rTi) est dans 
D{0,Roo) = -D(0, 0.005). Soit P2 de coordonnées T2, 92 tel que 52 est proche de 
qi, en ce sens que 92 € Z)(0,0.01). Pour tout o; e Pi on a \ Jg^ < |ç2 — 
qi\ max|,|<o.oi(| Efc>i Ofe?''"^]). Or pour |g| < 0.01 on a | J2k>i «fc?''"^! < Efe>o(fc+ 

6 

099^ 



1)^10-"'= <ô^<7. Donc 



|A*oo(f2) - Moo(A)|oo < 7\q2 -qi\ et |/Xoo(-P2) - f^oo{Pl)\2 < 7y/g\q2 - qi\ 

Soit maintenant Pi un point de coordonnées ti, qi et q'i tel que q[ ~ exp(- 



-2iTT - 



pri ' 

est dans D(Q,Rq) = -D(0, 1 — i). Soit P2 de coordonnées T2, q'2 tel que est 

proche de q[, en ce sens que € D{0, 1 — -^). Pour tout o; S Pi on a | J^/' uj\ < 

|Ç2 -ç'i|maX|çl<i__L(lEfe>iafc9*'"^l)- Or pour |ç| < 1 - ^ on a \T,k>iO-kq''~'^\ < 
Efc>o(fc + l)'(l-è)'^96p4. Donc 
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|Mo(f2) - Mo(Pi)U < 96/|g^ - q[\ et |/io(P2) - Mo(Pi)|2 < 96/V5k2 - «il- 
Ainsi la perte de précision occasionnée par l'application de Jacobi est O(logp). 

Lemme 3 (Majoration des intégrales de Jacobi) . — Pour tout premier p on 
pose Roo = 0.005 et Ro = 1 — ^ et on recouvre Xo{p) par les deux disques analytiques 

Doc et Dq centrés en chacune des deux pointes oo et et de rayons respectifs Roo 
et Ro. Donc D^ = D{0,R^) = {q^, \q^\ < 0.005} et Dq = D{0,Ro) = {qo, \qo\ < 

Sur chacun de ces deux disques, l'intégration de Jacobi définit une application 

à valeur dans le dual de l'espace des formes holomorphes TC^{Xo{p)). On munit 
H^^Xoip)) de la base Vi constituée des formes primitives, propres et normalisées 
et on note abusivement son dual, que l'on munit de la norme L°° associée à la 
base canonique (duale de Vi). 

L'applicaMon d'intégration de Jacobi est alors Lipschitzienne sur chacun des deux 
disques ( et même sur leurs voisinages D(0,0.01) et I?(0, l- ^)) et son coefficient 
de dilatation y est majoré par 7 sur le premier et par 96p'* sur le second. 

Soit e = {ek)i<k<g G {0,cxd}^. On note D^^ le produit 

D, = D{Q,R,,) X ■ ■■ X D{Q,R,^). 
Les De recouvrent le produit X(C)^. L'application produit 

associe au (/-uplet (çi, . . . , qg) la somme des ^i^^{q^^_k) pour 1 < k < g. 

On note {zuj)ui la base duale de la base canonique de C^^. On note jToi.e le 
déterminant jacobien de l'application /z^ en {qei,i, ■ ■ ■ , Q€g,g) '■ 



J-DiAlei,!, ■ ■ ■,Qeg,g) = " ^ il^iA' ■ ■ • ' = 



dqt 

Ce déterminant est une série entière de g variables dont les coefficients se ma- 
jorent aisément à partir d'une majoration des coefficients des formes lo. Mais il peut 
parfaitement s'annuler. U instabilité à'\m g-uplct dans peut se définir comme 
l'opposé du logarithme du module de ce déterminant jacobien. 

2.7. Jacobiens et wronskiens. — L'espace ti^ des différentielles holomorphes 
est muni de la base î?i et de la norme L°° associée. Si F = X^t^ePi cette 
norme est notée \F\ = \f \ = max^^ \fuj\- Puisqu'on dispose de deux disques analy- 
tiques sur X{€,) centrés en chacune des deux pointes, il est naturel d'introduire une 
norme S^{F) = max^^^sjo^^^) et de même Sq{F) = maXg^gS(o,H„) 

On introduit aussi les variantes <S(x)(i^) = maXg^g5(o,i) df~ et de même <So(F) = 

maXç^6D(o,i) -ifo ■ 

Puisque toutes les normes sont équivalentes, les quotients So{F)/\F\, Sao{F) /\F\, 
Sq{F)/\F\, Sao{F)/\F\ pour F ^ Q sont majorés et minorés par des bornes 
indépendantes de F. On veut montrer que le logarithme de ces bornes est polyno- 
mial en p. C'est évident pour la borne supérieure parce que les formes de la base 
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Vi s'écrivent uj = f{qx)dqoo = ±/((?o)c?<?o avec / d'ordre de grandeur^^^ (1)3). Il 
suffit d'appliquer le lemme 17 de majoration du reste. 

Pour contrôler la borne inférieure, il suffit de trouver g points qi, . . ., qg de 
£'(0, Roo) tels que le module du jacobien Jt>i,((x,...,(x){Q1j ■ ■ ■ ^Qg) ait un logarithme 
borné inférieurement par un polynôme — .g'^ en le genre g. 

Une stratégie possible est de chercher d'abord un q tel que le wronskien de Vi 
en q ne soit pas trop petit et de chercher ensuite des qi,. ■ -, qg dans un voisinage de 
ce q. 

On utilise le classique 

Lemme 4 (wronskien et jacobien). — Soient g >2 un entier naturel et fi{q), 
f2{q), ■ ■ -, fg{q) des séries de Laurent à coefficients complexes. On appelle wronskien 
associé à i = {fi, . . . , fg) le déterminant 



Wf{q) 



m) 



fgil) 



On se donne g indéterminées qi, q2, . ■ ., qg et on appelle jacobien associé à f le 
déterminant 



Jf 



fgin) 
fg{Q2) 

fgiQg) 



On note D = (9' ~ 9fe) discriminant réduit. 

Si les fk sont des séries entières alors le jacobien Jf est dans l'anneau 
€[[^1, . . . , Qg]] et il est divisible par le discriminant réduit D dans cet anneau. Le 
quotient Jt/D est alors congru à Wt{Q) modulo l'idéal maximal de C[[gi, . . . ,qg\[. 



La démonstration peut se faire par récurrence sur g. 



□ 



Lemme 5 (Majoration du wronskien). — Si pour 1 < l < g les fi{q) sont des 
séries entières d'ordre de grandeur {A, k) avec A > 1 et k > 1 alors le wronskien 
Wt{q) est une série entière d'ordre de grandeur 



{g\AS2 ê ^gk. 



1 



Cela résulte du lemme 15. 



□ 



Lemme 6 (Majoration du jacobien). — Si pour 1 < l < g les fi{q) sont des 
séries entières d'ordre de grandeur {A,k) avec A > 1 et k > 1 alors le jacobien 
Jt{qi, . . . ,qg) est une série entière des g variables qi, . . ., qg d'ordre de grandeur 



{g\As,{k,k,...,k)). 

Selon le lemme 4, le wronskien nous renseigne sur le terme dominant du jacobien 
au voisinage de gi = 92 = • • • = 9g = 0. 



f^^La définition de l'ordre de grandeur d'une série est donnée au paragraphe A.l. 
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Définition 1 (Wronskien et Jacobien des formes différentielles) 

Pour chaque forme différentielle primitive, propre et normalisée co GVi on rap- 
pelle que LO = ^^dq où f est une forme modulaire primitive, propre et normalisée 

de poids 2. On forme la famille f des g séries entières ainsi obtenues. Pour 
tout q G -D(0, 1); le wronskien de i en q est noté Woo{q) et le jacobien correspondant 
en q = (çi, . . . , Çg) est noté Joo(q)- 

Posons m = Le produit Woo{q){dq)"' est une forme holomorphe de degré 

m sur X{C). Donc elle a au plus 2{g ~ l)m zéros en comptant les multiplicités. En 
particulier, la pointe à l'infini q^o = est un zéro de multiplicité v < 2{g — l)m de 
W^{q). 

La série entière Woc (<z)^ a des cocfRcients entiers rationnels et sa valuation 2v est 
au plus 4((? — l)m = 2g{g^ — 1). On suppose > 2. Le lemme 5 de majoration du 
wronskien montre que l'ordre de grandeur de Woo{q)^ est {exp{g^), g^). Le lemme 
17 de majoration du reste permet d'écrire W^oiq)"^ = wq^^ -\- R2v+i{q) avec w entier 
naturel non nul et iÎ2D+i majoré en module par |gp''+^(2w+2)3 pour q G D{0, Roo)- 
Donc Woo{q)^ est minoré en module par ig^^^^^"^) si |g| < {2v + 2)~^ . 

Lemme 7 (Minoration du wronskien). — Il existe une constante effective pos- 
itive cq telle que pour tout premier p tel que le genre g de Xq{p) est au moins 2, il 
existe un q G Z>(0, 10"^'') tel que log|Woo(ç)| > —g'"^- Un tel q se calcule en temps 
polynomial en p et la précision absolue requise. 

Soit q G D{0, Roo) tel que Woo{q) 7^ et posons q = {q, . . . ,q) etx = (xi, . . . , Xg). 
Le jacobien j7'oo(q + {Ig — |lq|l)) est^^^ une série entière en les g variables xi, 
. . ., Xg. C'est en fait le jacobien des {fk{q + a^(l — \q\))/{q + — \q\)))i<k<g vues 
comme fonctions de x. Le lemme 4 appliqué à ces dernières donne le terme principal 
de cette série en x = Og : 

Joo(q + x^(lg - ||q||)) = Woo(q)(l - kl)^^ n(^' " Xk) + R,(,-i) , ^fx). 

k<l ^ 

Le jacobien j7oo(q) est une série en q d'ordre de grandeur {gl, (3, 3, . . . , 3)) d'après 
le lemme 6. Si g e D{0,Roç,), la série recentrée j7oo(q + x*(lg — ||q||)) est une série 

en X d'ordre de grandeur (cxp(g'-'), (4, 4, . . . , 4)) d'après le lemme 16 de recentrage. 
Pour X dans P{Og, Roo), le lemme 17 de majoration du reste donne alors pour g >2 

R a(a~i) I 1 (x) 

On pose s = |x|oo et on suppose que x = (|, y , . . . , '"^~g^^ , s) et s < iîoo- Alors 



Woo{<l){l-\q\)'^l[{xi-Xk) 
k<l 

> |Woo(q)| ( — ^ ) 
On choisissant un q donné par le lemme 7 et en prenant s = exp(— 3^) on prouve 

le 



< exp(5'^)|x 



+1 



> iw^ion{!^) 



f^^La notation ★ est introduite au paragraphe A.l. 
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Lemme 8 (Minoration du jacobien) . — Il existe une constante effective 
cr telle que pour tout premier p tel que le genre g de Xo{p) est au moins 
2, il existe r = (ri,...,rg) G D{0,R^y tel que log|Joo(r)| > —g"'^- Un tel 
r se calcule en temps polynomial en p et la précision absolue requise. Les 
cinq normes sur définies pour F = J2ueT>i f<^^ P^''' \^\ ~ max^^ 1/^1» 
Soo{F) = max gg(o,floo) d£: ' ^o{F) = max g_D(o p„) -£- , S^{F) = 



max„ 



= e-D(o, 



j— , et So{F) = maXgpg5(Q 1-, ^ sont dans des rapports 
exp(p^), ou, si l'on préfère, il existe une constante positive C3 telle que pour 
tout p et pour tout F ^ les différences entre log|F|, logiSo(-F), Iog5oo(-F), 
log<So(-F) etlogSaoiF), sont bornées en valeur absolue parp^^. 

2.8. Jacobiens et wronskiens quadratiques. — Des résultats semblables sont 
vrais pour le jacobien et le wronskien associés à la base de l'espace de formes 
différentielles quadratiques W^(A2). 

Définition 2 (Jacobien et wronskien des formes quadratiques) 

Pour chaque forme différentielle quadratique primitive, propre et normalisée (j) G 
on rappelle que = blâl'éÉ— /j ^gi yjjg forme modulaire primitive, propre 
et normalisée de poids 4. On forme la famille h des 52 = 3gf — 1 + f 2 + î^s séries 
entières ainsi obtenues. Pour tout q G -D(0, 1), le wronskien de h en q est noté 
W2,oo il) et le jacobien correspondant en q = (91 , • • • , ) ^oté 1/2,00 (q) • 

Posons TO2 = 92(92+3) ^ prodnit W2, 00 (q) q~ (dq)"^^ est une forme de degré 
m2 sur X{C), holomorphe en dehors de A2. Plus précisément elle appartient à 
H"*^ (52 A2). Donc elle a au plus 2{g — 1)to2 + 52(2 + 1^2 + l'a) zéros en comptant les 
multiplicités. En particulier, la pointe à l'infini q^ = est un zéro de multiplicité 
V2 < 2{g - l)m2 + ,92(2 + 1^2 + v-i) de W^2,oo('7)- Donc V2 < 9g{g + 1)^. 

La série entière W2^oo{q)'^ a des coefficients entiers rationnels et sa valuation 2i>2 
est au plus 18g(g + l)^. On suppose g >2. Gommes les sont d'ordre de grandeur 
(1,4), le lemme 5 de majoration du wronskien montre que l'ordre de grandeur de 
W2.oo{q)'^ est (exp(5'^), g'-'). Le lemme 17 de majoration du reste permet d'écrire 
^^2,00(9)^ = wq'^'"^ + R2v2+iiq) avec w entier naturel non nul et R2V2+1 majoré en 
module par \q\'^'"^~^^ (2v2 + 2)^ pour q G D{0,Rao)- Donc ^2,00(9)^ est minoré en 
module par 

iqisg{g+ir si -log|g| > 5^ et 5 > 2. 

Lemme 9 (Minoration du wronskien quadratique). — // existe une cons- 
tante effective cg telle que pour tout premier p tel que le genre g de Xçi(p) est au 
moins 2, il existe un q G D(0, 10^^°) tel que log |W2,oo(q')| > —g'^^- Un tel q se 
calcule en temps polynomial en p et la précision absolue requise. 

On obtient de même l'analogue du lemme 8 pour le jacobien des formes quadra- 
tiques. 

Soit q G D{0,Roc) tel que ^^2,00(9) 7^ et posons q = {q,...,q) et x = 
{xi, . . . , Xg^). Le jacobien J72,cc (q+x*(lg2 — ||q||)) est une série entière en les 92 vari- 
ables xi, . . ., Xg^. C'est en fait le jacobien des {hk{q+x{l—\q\)) / {q+x{l — \q\)))i<k<g2 
vues comme fonctions de x. Le lemme 4 appliqué à ces dernières donne le terme 
principal de cette série en x = Og^ : 



92(92"-'-) T — r 

J2,oo(q + x) = VF2,oo(q)(i - \q\) = [[(xi - a^fc) + -^ 92(32-1) , i(x). 

k<l 



JACOBIENS, JACOBIENNES ET STABILITÉ NUMÉRIQUE 



13 



Le jacobicn J2,oo(q) est une série en q d'ordre de grandeur (52!, (4, 4, . . . , 4)) 
d'après le lemme 6. Si g G D(0,iîoo)) la série recentrée »72,oo(q + x* (Ig^ — ||q||)) 
est une série en x d'ordre de grandeur (exp((7^), (5, 5, ... , 5)) d'après le lemme 16 
de recentrage. Pour x dans P{Og,Roo), le lemme 17 de majoration du reste donne 
alors pour g >2 



82(92-1) I ^ 

-^ 82(32-1) I ^(x) < exp(5 )|x|oo ^ 



On pose s = |x|(x> et on suppose que x = |f , • • • 



(92-1)^ 



, s) et s < iîoo . Alors 



k<l 



> \W2,oo{^)\ 



> \W2,M\ 



5(1 -M 



92 



0.995s 
92 



82(82-1) 



82 (82-1) 



On choisissant un q donné par le lemme 9 et en prenant s = exp(— g^) on prouve 

le 

Lemme 10 (Minoration du jacobien quadratique). — E existe une cons- 
tante effective cg telle que pour tout premier p tel que le genre g de X(){p) est au 
moins 2, il existe r = {n, . . . jrg^) G £'(0, iîoo)^^ tel que log 1 1/2,00 (r)| > —g'^^- 
Un tel r se calcule en temps polynomial en p et la précision absolue requise. Les 
cinq normes sur W^(A2) définies pour H = '^^^Dz ^'('^ P'^^ \^\ ~ max^i 



Soo{H) = max, 
°iax,^eS(o,i) 



gcx=e-D(0,floo) 



et log Sf) [H) sont bornées en valeur absolue par p'^^ . 



"^9'Xgo6£'(0,flo) 

qgH 



go g 



Soo{H) 



, sont dans des rapports 



, et So{H) = max^^g5(o_i) 
exp(p^) ou, si l'on préfère, il existe une constante positive 04 telle que pour tout p 
et pour tout H ^ les différences entre log \H\, log»So(iî), \ogSoo{H), log»Soo(-ff), 



2.9. Stabilité. — On suppose que g > 2. On suppose choisie une origine o et on 
note S^X la g-ième puissance symétrique de X et : S^X C^i/£ l'apphcation 
d'intégration de Jacobi. Soit e = {€k)i<k<g € {0,oo}^ et r = (n, ■ ■ ■ ,rg) G = 
D{0, iîe J X • • • X D(0, Reg) et notons /? e l'image de r par fj,^. On rappelle que 
l'instabilité de r est l'oppose du logarithme du module du déterminant jacobien 



Jvi, 



. On suppose cette instabilité finie et on la note A. Donc la 



restriction de ji^ à un voisinage de r est injective. L'image d'un tel voisinage est un 
voisinage de p et la corestriction de fjp à un voisinage assez petit de p est injective. 

On veut étudier quantitativement cette situation. On montre d'abord que les points 
voisins d'un point stable sont assez stables. Ensuite on montre que l'image par 
l'application de Jacobi d'une boule centrée en un point stable, contient une boule 
pas trop petite. 

D'après le lemme 6, le jacobien Jtix,^ est une série en les Çe^.fe d'ordre de grandeur 
(5!, (3, . . . , 3)). Posons X = (xi, . . . , Xg). La série recentrée J2?i,e(i" + x* (Ig — ||r||)) 
est d'ordre de grandeur (exp(gi^), (4, ... ,4)). Pour |x|oo < \, le reste d'ordre 1 de 
cette série est majoré par exp(g'^)|x|oo. Donc si |x|oo < exp(— — A) le point 
r + (Ig — ||r||) ★x est stable d'instabilité majorée par 1 + A. Si y = (Ig — ||rl|) *x 
alors |x|oo < p|y|oo et |y|oo < |y|2 donc si |y|2 < exp(-g'^ - A) le point r + y est 
stable d'instabilité majorée par 1 + A. 
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On compare maintenant l'application do Jacobi j.!, et sa différentielle D/Ze en 
r. La matrice de D/i^ dans les bases {dqf^)k et {zi^)i^, la base dualc de la base 
canonique de C^^, n'est autre que (^^(^fe)) • EjUc a des coefficients majorés 

par 3l2p'^. Donc pour y 7^ Og le quotient des normes L°° associées aux deux bases 
susmentionnées ^^1^' ^ ^^°° est majoré par g^. Puisque le déterminant de Dfj,^^ est 

minoré en module par exp(— A) on montre de même que le quotient ^^i^,^^"" est 

|y|oo 

minoré par exp(— A)5f~^^ > exp(— A — Og"^) donc 



(2) exp(-A - Og^)\y\2 < \Dfi,y\2 < g°\y\2. 

Soient u une forme dans Vi et k un entier entre 1 et g. Le développement de 
uj sur = D{0,B.e^) s'écrit uj = i^^^^^dg^^ pour une forme primitive, propre 
et normalisée / de poids 2. La série f{qek)/(lek ^ pour ordre de grandeur (1,3). 
Sa recentrée en a pour ordre de grandeur {g^,4). Si ç^j, = rfe + — |rfe|), 
avec < ^, le reste d'ordre 1 de la série recentrée est majoré par g^\xk\- Donc 
jr.+^.d-r.) ^ _ ^^^^ _ ^^^f^ ^ ^o^^ _ rk)\xk\^. Ainsi le reste d'ordre 2 



iÎ2Me((l9-r)*x) =/i,(r+(lg- |lr||)*x)-/Xe(r)-7^Aie((lg- ||r||)*x) 

est majoré en norme L°° par g(^|x|^ et en norme par £f^|x|| < ff^lyll avec 
y = (Ig - ||r||)*x. Ainsi 



(3) \R2f^e{y)\2 = lA*.(r + y) - Me(r) - Df,,{y)\2 < g'^WH 

Les équations (2) et (3) permettent de s'assurer que la partie principale D/i^y est 
deux fois plus grande que le reste R2n{y) : si — log|y|2 > Og"^ + A alors \D^fy\2 > 
2|iÎ2M(y)|2 donc /x,(r + y) - /Ze(r) > i|D/i,y|2 > exp(-A - Og^)\y\2. 

Ainsi pour tout 6 > g''^ , l'A boule centrée en p et de rayon exp(— Og^ — 2A — 6') 
pour la norme est formée de points stables (ayant un antécédent unique par 
l'application de Jacobi) et elle est contenue dans l'image de la boule centrée en 
r et de rayon exp(— A — 9). 

Lemme 11 (Stabilité du problème inverse de Jacobi) 

Il existe une eonst,ante effective positive Cio telle que l'énoneé suivant soit vrai : 
Soit e = {tk)i<k<g dans {0, oo}'' et v = {n, . . . ,rg) € = £)(0, iî^J x • • • x 
D{0,Re^) et notons p l'image de r par l'application d'intégration de Jacobi p^- 
On suppose que l'instabilité X de r est finie. Alors pour tout 6 > g'^^° la houle 
centrée en r et de rayon cxp(— A — 6) est form,ée de points stables et son image par 
l'application de Jacobi contient la boule centrée en p et de rayon exp{—g'^"' — 2X — 9). 
Donc si r' G S^X s'envoie sur p' par l'application de Jacobi et si \p — p'\2 < 
exp(-25'=i'' - 2A) alors |r' - r|2 < exp(5f^io + A)|p - p'\2. 



2.10. Quelques sous ensembles discrets de la jacobienne. — On se place 
à nouveau dans la situation du paragraphe précédent. On suppose que g > 2. On 

suppose choisie une origine a et on note S^X la g-ième puissance symétrique de X 
et p^ : S^X C^'-fC l'application d'intégration de Jacobi. Soit e ~ {ek)i<k<g € 
{0, oo}3 et r = (n , . . . , Tg) e De = .0(0, Re^ ) x ■ ■ ■ x D{0, R^^ ) et notons p l'image 
de r par p^. On note A l'instabilité de r. 
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On identifie D -D^ à l'espace tangent en r et on note {Sk)i<k<g la base duale 
de {dqe^)i<k<g. 

L'espace vectoriel réel sous-jacent a pour base (^1,^2, ■ ■ ■ ,Sg, iôi, iÔ2, ■ ■ ■ , iôg). 
Pour k entre 1 et 5 on pose iôk = Sk+g- 

On choisit un réel x > logp, on pose T = cxp(— x) et on considère les 2g 
petits accroissements . . ., (32g de C^/£ définis par (3k = ^^{rk + T) — /i(rfe) = 
(/rfc'°"^^^)<-'625i = POur l<k<get (3g+k = Kvk + i'^)- ^Vk) = 

ilr^^^'^ ^)u>eT>r = + '^ih) - /"e(r) pour 1 < fc < 5f. 

Si M est un entier positif, on note A{r,x,M) le sous ensemble de C^/C formé 
des combinaisons des /3k à coefficients entiers dans [—M, M]. On se demande dans 
quelle mesure les éléments du tore C^/C peuvent être approchés par l'ensemble 
A{r,x,M). 

On doit premièrement minorer le déterminant f^^'^^^f^" où (6^)^ est la base 
canonique de l'espace réel sous-jacent à = En particulier (efe)i<fc<g 

est la base canonique de et Ck+g = ick- On note D/j^ la différentielle de en r 

7iA---A72 
ei A---Ae2 



Pour 1 < A; < 2^ on note 7/5 = Dfj,J5k)- Le déterminant ''^^'"^''^s est le carré du 
module du déterminant jacobien. 

Or T7fe est une bonne approximation de (3k- Plus précisément, on a vu au para- 
graphe précédent que si — logT = x > Og"^ + A alors T|7fc|2 > 2|/3fe — Tji;\2 donc 
il7fe|2 < Wk\2 < ÇllfcU < d'après l'inégalité (2). 

L'inégalité (3) implique quant à elle que \(3k — "^JkU < ff'^'ï"^- 

Comme le déterminant est multilinéaire, on majore le module de la différence 
TT^^ - T^'^irX^ par (maxi<,<2, l/J.b)'^"' niaxi<fc<2, |/3, - TjkU- 

Ainsi la différence ^l^UlM^a _ j2g iiA-A72, majorée par 22350(23-1)^0x23+1. 
Cette différence est plus petite que la moitié de T^^'l l'X-Tell l = exp(-2A)T2s 
pourvu que x > + 2A. Alors 

— — ^ > 7; exp -2A - 2gx). 

ei A • • • A e2g 2 

On a donc une minoration du déterminant de la matrice de passage de la base 
(6fc)i<fc<2g à la base {(3k)i<k<2g- Les colonnes de cette matrice sont les (3^ et elles 
sont majorés en norme L°° par g'~'T. Donc les colonnes de la matrice inverse sont 
majorées en norme L°° par 2(2^ - 1)! exp(2A + 2gx)g^^'^^~'^^T'^^~'^. 

D'après le lemme 2 il existe une constante positive C2 et une base du réseau C 
des périodes de coordonnées majorées en valeur absolue par exp{p^^) dans la base 
canonique de R^^. Donc les coordonnées de ces périodes dans la base {Pk)i<k<2g 
sont bornées par 2{2g - 1)! exp(p'== + 2A + 2gx)g^^'^'^~'^^T^'^~'^ < exp(2A + x + 5^)- 

Les coordonnées de tous les points du parallélogramme fondamental associé à 
cette base sont également bornés de la sorte. Enfin, en remplaçant de telles coor- 
données par leurs valeurs entières on commet une erreur bornée en norme L'^ dans 
la base canonique par g^T. 

Lemme 12 (Bonne répartition). — E existe une constante C5 effective positive 
telle que : si e G {0, 00}^ ei r G a une instabilité finie \ et si x est un réel plus 

grand que g'^'^ + 2A et si M est le plus petit entier plus grand que exp(g'^^ + 2A + x), 
alors tout point p de = IR.29 ggi à distance < exp(— x + g'^^) d'un point de 
A{r, X, M) et les coefficients entiers dans [—M, M] de ce point peuvent être calculés 
en temps polynomial en g, log(max(l, IpIoo)) et x- 
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3. Complexité des opérations dans la jacobienne 

Dans le paragraphe 3.1 on étudie la complexité de l'algorithme qui, étant donné 
une famille de points, trouve une forme s'annulant en ces points, puis recherche les 
autres zéros de cette forme. On en déduit dans le paragraphe 3.2 que les opérations 
élémentaires dans Jq (p) sont stables et se font en temps déterministe polynomial en 
p et la précision requise. 

Au paragraphe 3.3, un semblable résultat est établi pour le problème inverse de 
Jacobi. 

3.1. Une opération élémentaire dans la jacobienne : la dualisation. — 

Un point P de X est spécifié par une valeur de r G H* ou de qo G -D(0, Rq) ou de 
Qoo € D{0,Roo)- Connaître un tel point, c'est disposer d'un algorithme qui calcule 
par exemple le Qoo correspondant, en temps polynomial en g et la précision absolue 
requise. 

On suppose g > 2 et on se donne 3g — 4 points distincts Ri, R2, . . ., Rsg-A sur 
X. On suppose que ces points ne sont ni des pointes, ni des points elliptiques. On 

pose e = 3.g — 2 + + 1^3 = 52 — 1 et on note Pi , .... Pp la famille de points formée 
des 3(7 — 4 points que l'on vient de se donner et des 2 + 1^2 + 1^3 points elliptiques ou 
paraboliques. 

On suppose que les r premiers points Pi, . . ., P,. sont dans Dqo et sont donc 
donnes par des goo,/c € D(0,Roc) pour 1 < k < r. Les autres points sont supposés 
appartenir à Dq et sont donc donnés par des go,fe G -^(0, Rq) pour r + 1 < k < e. 
On suppose que Pi est la pointe 00 et Pg la pointe 0. 

On se donne un 77 > et pour tout k entre 1 et r on considère le disque D^ k = 
{looi \<loo — 9oo,fe| < de rayon -q centre en Çoo.fc- De même, pour tout k entre r + 1 
et e on note Do,fe = {<Z0) \lo — 1o,k\ < v} le disque de rayon 77 centré en go,fc- On 
suppose que les images de tous ces disques dans X{C) sont deux à deux disjointes. 

Définition 3 (Résolution d'une famille de points). — Dans cette situation, 
on dit que {{Pk)i<k<Sg-2+v:i+v3-i'i',ri) est une résolution pour la famille de points 

{Rk)l<k<3g-i- 

On note P = Pi + • • • + P(. le diviseur somme des Pj. et on observe que P + 
A2) C H^(A2) est un sous-espace vectoriel non nul de Ti? . 

On représente une forme F dans H^(A2) par ses coordonnées dans la base D2 soit 
F = X^^ePa •^<^'^' cherche les équations linéaires qui définissent le sous espace 
Ti?{—P + A2) dans cette base. Il y a une équation pour chaque Pfe. 

L'équation correspondant à la pointe à l'infini Pi, est 

°^ {dq\^ ^ X]fe>i '^<t>-k<llL est le développement en l'infini de la forme primitive, 
propre et normalisée de poids 4 associée à 4>. 

L'équation correspondant à la pointe en zéro Pg, est 

(dqo)'^ ~ X]fe>i^0,fc9o le développement en zéro de la forme primitive, pro- 
pre et normalisée de poids 4 associée à (p. Noter que 6^,1 = ±a^,i car les formes 
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primitives, propres et normalisées sont en particulier valeurs propres de l'involution 
d'Atkin-Lehner. 

L'équation correpondant à Pfc pour 2 < fc < r est 



(dq\^ (Pfc) = Xlj>i '^<P,illx \ calcule en temps polynomial en p et la précision 
requise. 

L'équation correpondant à Pk pour r + l<fc<e— lest 



On rassemble toutes ces équation dans une matrice M.p dont le noyau décrit 
7{^{—P + A2) dans la base V2. Les coefïicients de cette matrice sont < p'^ . On 
calcule une valeur approchée A4'p de A^p à coefficients dans Q{i). On note m > 2 
la précision de cette approximation. Donc M'p — Mp a, des coefficients inférieurs à 
exp(— m) en module. Il existe une telle A4'p de taille polynomiale en log^> et m et 
on la calcule en temps polynomial en p et m. L'algorithme d'Hermite [11], ou mieux 
l'algorithme LLL [2, Theorem 2.6.2], produisent une base à coefficients entiers du 
noyau de A4'p de taille polynomiale en p et m en temps polynomial en p et m. On 
choisit un vecteur de cette base. Si M = a + ib est le plus grand des coefficients 
de ce vecteur, on normalise en divisant le vecteur par le maximum de \a\ et |6|. On 
obtient ainsi un /' = {f'^)(j,eTi2 dans le noyau de M'p à coefficients dans Q{i) et 
F' = f'^(j) la forme associée telle que \F\ = max^ |/^| soit compris entre 1 et \/2. 
D'après le lemme 10, les normes So{F') et S^{F') sont majorées par \/2 exp(p'^'' ) 
et minorées par cxp(— p'^*) pour imc certaine constante positive C4. 

Le produit Mpf = {Mp — M'p)f' est un vecteur de coefficients majorés en 
module par 52 exp(— m)v^. Ces coefficients sont, dans l'ordre, la valeur de i^J^^yA 
en la pointe infini, les valeurs de i^J^^y en les Pk pour 2 < fc < r, les valeurs de 

i^J'qgyi en les Pk pour r + l<fc<e — 1, la valeur de (f^yi en la pointe nulle Pe- 
On peut maintenant appliquer le lemme 21 de stabilité des zéros. 
On s'intéresse d'abord aux P^ pour k < r. Ils sont contenus dans £'(0, iîoo)- 
Puisque < V^, la série entière çJî^^y est d'ordre de grandeur (172 V^, 4). Avec 
les notations du lemme 21 on a donc A = gi\f2^ n = 4, 1 — |c| > 0.995 et — loge = 

(0.995)^'* /m" 

m/2 pourvu que m > 01og(g(). On demande que — — ^ ^ soit plus grand que 
Ci4(4^(l — log0.995)+log5(2-\/2) et quep'^* ce qui est assuré si m > p^. Alors chaque 
Pfc pour A; < r est à distance < exp(— A/m/2) d'un zéro de F' . 

On s'intéresse maintenant aux P^ pour k > r. Ils sont contenus dans D{0, Rq) 
avec Ro = 1 — -■ Puisque < \/2, la série entière est d'ordre de grandeur 

(5'2\/2, 4). Avec les notations du lemme 21 on a donc A = gi\f^, = 4, 1 — jcj > | 

et — loge = m/2 pourvu que m > 0\og{g). On demande que soit plus grand 

que ci4(4^(l + logp) +log52\/2) et que p'^^ ce qui est assuré si m > p^. Alors chaque 
Pfc pour k > r est à distance < exp(---\/m/2) d'un zéro de F' . 

On s'assure enfin que exp(— •\/m/2) est plus petit que la résolution 77 de sorte 
que les e zéros ainsi trouvés sont distincts. On cherche alors les g autres zéros 






(-ffc) = J2i>i ^<l>,i1o k ^® calcule en temps polynomial en p et la précision 
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de F' avec l'algorithme de quadrichotomie de Weyl (voir le lemme 22). Puisque 
l'application d'intégration de Jacobi est Lipschitzienne avec coefficients de dilatation 
polynomiaux en p on obtient le 

Lemme 13 (Dualisation). — Il existe un algorithme déterministe qui, étant don- 
nés un nombre premier p et 3g — 4: points distincts Ri, . . ., R^g-i sur Xçj{p) et 
une résolution {{Pm)i<m<g-2T f'^'l)- calcule g points Qi, Qg dont la somme est 
linéairement équivalente à deux fois le diviseur canonique moins la somme des points 
initiaux : 

Ql + • • • + Qg - 2/C - (iîi + • • • + iÎ3g-4). 

Cet algorithme est polynomial en p, le logarithme — log de la résolution et la 
précision requise. 

Remarque importante : Notons que dans cet énoncé, l'erreur qui est majorée 
se mesure dans l'image C^^ de l'application de Jacobi. 

On peut s'affranchir de la condition sur la resolution r/. En effet, si les points 
Rk ne sont pas tous disctincts, ou plus généralement si leur résolution est jugée 
trop faible, c'est à dire si r] est trop petit, on peut perturber légèrement ces points 
pour obtenir la résoluton souhaitée. Comme l'application d'intégration de Jacobi 
est Lipschitzienne de dilatation polynomiale en p, on peut adapter la perturbation 
à la précision finale souhaitée. On obtient ainsi le 

Lemme 14 (Dualisation). — Il existe un algorithme déterministe qui, étant don- 
nés un nombre premier p et 3g — A points Ri, . . ., Rsg-i sur Xq{p), calcule g points 
Qi, Qg dont la somme est linéairement équivalente à deux fois le diviseur 
canonique moins la somme des points initiaux : 

Ql + • • • + Qg ~ 2/C - (iîi + • • • + iÎ3s-4). 

Cet algorithme est polynomial en p et la précision requise. 

3.2. Addition et soustraction dans la jacobienne. — On suppose que le 
genre de Xq{p) est au moins 4. On choisit un diviseur effectif origine de degré g 
noté O = Oi + • • • + Og . On choisit aussi un diviseur effectif auxiliaire de degré g — 4 
noté N . La dimension de 7i^(A2) est 92 = 3,g — 1 + + '^a- Le plus naturel est 
de choisir un point o origine de l'application d'intégration de Jacobi et de poser 
O = go et N = {g — 4)o. Un élément de Pic°(X) est donné comme classe d'un 
diviseur Q — O oii Q est effectif de degré g. Soit R un autre diviseur effectif de 
degré g. Pour ajouter la classe de Q — O et celle de iî — O on applique le lemme 14 
au diviseur Q + iî + TV qui est effectif de degré 3g — 4. On obtient un diviseur T 
effectif de degré g tel que T 21C — Q — R — N . On applique à nouveau le lemme 
14 au diviseur T + O + A'^ et on obtient un diviseur U effectif de degré g tel que 
i7 + O ~ Q + iî. Donc [/ - O est bien la somme de Q - O et iî - O. 

Pour calculer l'opposé de Q — O on applique le lemme 14 au diviseur 20 -\- N ce 
qui produit un diviseur S effectif de degré g équivalent &21C — N— 20. On applique 
à nouveau le lemme 14 au diviseur S + Q + et on trouve un diviseur V effectif 
de degré g tel que V — O ^ —{Q — 0). 

Théorème 1 (Arithmétique dans la jacobienne). — Les opération d'addi- 
tion et de soustraction dans la jacobienne de Xq {p) se font en temps déterministe 
polynomial en p et la précision requise. 
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3.3. Le problème inverse de Jacobi. — Étant donné p = {pui)u) € on 

cherche un diviseur effectif de degré g sur X qui s'envoie sur p modulo le réseau C 
des périodes par l'application de Jacobi. 

Le lemme 12 de bonne répartition, que l'on applique au r fourni par le lemme 8 de 
minoration du jacobicn, et à un réel x assez grand, fournit un ensemble ^(r, M) = 

{t\l3\-\ \-t2gP2g avec —M <tk< M} où les Pk sont les images par l'application de 

Jacobi de diviseurs Qk connus, différences de deux points proches. Les coordonnées 
réelles de p dans la base formée des pk sont calculées en inversant une matrice 
dont le déterminant n'est pas trop petit, puis tronquées à l'entier inférieur ou égal 
le plus proche. Les entiers tk ainsi obtenus, on calcule dans Jo(p) la combinaison 
J2i<i<2g ik^k avec la méthode du théorème 1 d'arithmétique dans la jacobienne et 
en utilisant l'exponentiation rapide. Le nombre d'opérations élémentaires dans la 
jacobienne est donc polynomial en logM, donc aussi la perte de précision. 



Théorème 2 (Problème inverse de Jacobi). — Étant donné un entier pre- 
mier p, on note g le genre de Xq{p) et Vi l'ensemble des formes différentielles 
primitives, propres et normalisées et C le réseau des périodes de Xq{p). 

Il existe un algorithme déterministe qui pour p € et pour O = 0i-\-- ■ ■-\r0g une 
origine de l'application d'intégration de Jacobi p^ : S^Xq{p) C'^'^/C et pour k un 
entier positif, calcule g points Pi, . . ., Pg de Xo{p) tels que \p^{Pi, . . . , Pg) — p\^ < 
exp(— fc), en temps polynomial enp, la précision k, et la taille log(max(l, IpIoo)) de 
P- 



On note que la dépendance en log(max(l, |/3|oc)) est sans importance car le reseau 
des périodes admet un parallélogramme fondamental de rayon exp{p'-') d'après le 
lemme 2. 

Notons encore que dans cet énoncé, l'errcnir qui est majorée; se; mesure dans 
l'image de l'application de Jacobi. Si l'on veut majorer l'erreur commise sur les 
points Pfe on note Qi, . . ., Qg des points tels que p'-'iQi, ■ ■ ■ , Qg) = p et soit A leur 
instabilité que l'on suppose finie, sans quoi la question serait dépourvue de sens. 

Le lemme 11 de stabilité du problème inverse de Jacobi, montre que la perte de 
précision dans le problème inverse de Jacobi est polynomiale en p et A. Il s'agit donc 
de contrôler A. 

Soit e G {0, oo}3 tel que Qk € D^,^ pour tout 1 < k < g et notons q^^^k la valeur 
de qc^. en Qk- L'instabilité A est l'inverse du logarithme du module du déterminant 



jacobien 



dq (<?£fc,fc) 



qui n'est pas une fontion algébrique. Il est donc naturel de 



réécrire ce jacobien comme le produit d'une quantité algébrique et de facteurs plus 
simples. 

Pour tout k entre 1 et g, on suppose que la valeur j{Qk) de l'invariant modulaire 
j en Qk vérifie j{Qk) ^ {0, 1728, 00}. 

On note que j{Qk) ne peut être proche de et 1728 en même temps. Si j{Qk) n'est 
pas proche de on pose jk = j — 1728 et sinon on pose jk = j- Le jacobien se réécrit 
alors comme produit du déterminant i/aig = ^^{Qek,k) et des . {Qek,k) pour 

J aJ,fc k 

l<k<g. 

Supposons que les Qk sont des points algébriques sur Q et que l'instabilité est 
finie. Alors Jaig est fini et non nul et on peut le minorer en fonction de la hauteur 
des j{Qk)- Les j^^g {lek,k) se minorent en utilisant le lemme 22 de stabilité des 
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zéros d'une série entière^^^ : si -^—{qe^,k) est petit alors q^^^k est proche d'un zéro 
de . Mais le jk choisi est alors plus petit encore. 

Appendice A 
Appendice sur les séries entières 

On a rassemblé dans cette section les définition et résultats relatifs à la local- 
isation et à la stabilité des zéros des séries entières qui sont nécessaires à notre 
travail. 

Le paragraphe A.l introduit quelques définitions et notations ainsi que des 
résultats élémentaires. Le paragraphe A. 2 énonce et démontre une forme simple 
et quantifiée du théorème de prolongement analytique. On introduit le polygone 
de Newton d'une série entière dans le paragraphe A. 3 et on montre comment il 
permet de localiser les zéros de cette série. Le paragraphe A. 5 prouve un résultat de 
stabilité des zéros d'une série entière et en déduit une majoration de la complexité 
de la localisation des zéros. 

A.l. Ordre de grandeur. — Soit 5 > 1 un entier. Si x = {x\, . . . ,Xg) G C^, 
on note |x|oo = max^ \xk\ la norme L°° de x. On note |x|i = \xk\ la norme 
de X et |x|2 = \/J2k l^^l^ norme L^. On note ||x|| le vecteur (|.ti|, . . . , \xg\). Si 
y = (2/1, . . . ,î/g) on note x*y le vecteur {xiUi, . . .,Xgyg). On note Og = (0, . . . , 0) G 
et Ig = e C9. Si X = {xi,...,Xg) e M», on dit que x > Og si et 

seulement si xj. > pour tout k. On dit que x > Og si et seulement si x^. > pour 
tout k. Notons P(x, r) = nfe=i D{xk,rk) C le polydisque de polycentre x et de 
polyrayon r. 

Une série entière / est donnée par / = /kx"* oii k parcourt N®. 

Définition 4 (Ordre de grandeur). — Soit A>1 un réel etn= (m, . . . , rig) G 
tel que n > Ig. On dit que f est d'ordre de grandeur {A, n) si pour tout k > Og 
on a 

|/k| < A(k + lg)" = A II {km + ir-. 

Si f est d'ordre de grandeur {A, a) et h d'ordre de grandeur {B, b) alors le produit 
p = fh est d'ordre de grandeur {AB, a + b + Ig). 

En effet p = J^mPm^"' ^^ec Pm = Ek+i=m /k'î'i- Il y a ni<„<g("^n + 1) = 
(m + Ig)"*"' termes dans cette dernière somme, et chacun est majoré en module par 
\hhi\ < A{k + lgTB{l + lg)b < AB{k + 1 + lg)'^+b. □ 

Lemme 15 (Dérivée). — Si f est une série entière d'une variable, d'ordre de 
grandeur {A, m), alors sa dérivée d-ième est d'ordre de grandeur (A2'''""' , to + 
d). 

Soit 5 > 1 un entier et / = X^k /k^'' une série entière de g variables d'ordre de 
grandeur (^,n). On en déduit pour tout z G -P(Og, 1) la majoration 



f^^ou un argument de compacité, plus général mais non effectif. 
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k>Op i<m<g km>0 

n\A 

< 



n\A n\A 

< 



(Ig- ||z||)°+i« - (l-|zU)9+Hi' 
Soit k = {kl, . . . ,kg). Posant pour tout m 



fcm + 1 + {Um + 1)1 



l'intégrale de Cauchy donne 



(2 

< An!k! 
= An!k! 



— f [ [ f{CiX2,...,Q ^^^^^^ 



(k + n + 2g)''+"+29 1 

(k + l,)k+l«(n + nrr>(l - \Zm\r-+''-+' 

(k + n + 2g)''+"+2p 



• (k + („ + l,)n+i., (1, - ||z||)"+>'+2. 

(4) < An!k!- 



(k + n + 2g)''+"+29 



■ (k + lg)^+^<' (n + - \z\^)'^a+\n\i+\kW 

Etant donné c = (ci, . . . , c^) S -P(Og, 1) il est alors naturel de considérer la série 
recentrée en c 



Fc{y) = Fc{yi,...,yg) = /(c + y*(lg - ||c||)) = /((c„ - \cm\))m) 

définie pour y G P{Og, 1). 

Reprenons l'inégalité (4) et déduisons 

Lemme 16 (Recentrage). — Soit g > 1 un entier, A > 1 un réel et n > Ig 
dans et soit f — X]k>Og f^'^^ ^'^^ série entière d'ordre de grandeur {A,n). Soit 
c G P{Og, 1) et notons Fc{y) = /(c + y * (Ig — ||c||)) la série recentrée de f en c. 
Alors Fc est une série entière d'ordre de grandeur (^c,n+ Ig) avec 



A, = n!Aexp(5 + |n|i)2«+l"lHlg - ||c||)-"-^ 
Pour tout entier positif u on note Ru{^) le reste de la série 



/(z)= /kz''+iî„(z). 

|k|i<u-l 

Soit 0<iî<letze P{Og, R) on a 
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\Ru{^)\ 



ikii>« 



< E E 1^^' 

|k|i>u |k|i>u 

k<(tt-l)lg k^(u~l)lg 



' (1- |z|oo)9+l"li 

avec n = |n|<^ et B = ^^^^â^. 

Lemme 17 (Majoration du reste). — Soit /(z) une série de g variables d'or- 
dre de grandeur {A,n). Soit u un entier positif ou nul et iîu(z) le reste d'ordre u. 
Soit R un réel strictement compris entre Q etl. Pour tout z dans P{Og,R) on a 

|iî„(z)|<i?(« + l)("+i)«|z|^ 



avec n = \n\^ et B = ^^jf^^^M, ■ 

En outre, siQ < k <1 est un réel et si 



,lQ{ngf 2(logK-logB) 

u > max -7T, -) 

^(logiî)2' logiî ' 



alors |-Rm(z)| < n pour 2, G P{Og,B} 



En effet, si w > alors {n + l)g^. < HiMM donc logS + (n+ l)c/log(l 

^<logKcaru>^ii2£^ 



u) + u\ogR< log B + < log KC&Tu> 2(i°g «-i°g B) _ □ 



A. 2. Prolongement analytique sur un disque. — Dans ce paragraphe on 

veut montrer qu'une série entière / d'une variable et d'ordre de grandeur (A, n) avec 
A> 1 et n > 1, majorée en module par un e > sur un petit disque D = D{c, r) 
inclus dans D{1),1) peut être agréablement majorée en module sur le gros disque 
D{0,^). 

On procède par prolongement à des disques de plus en plus gros. 
On introduit donc la 

Définition 5 (Fils d'un disque équilibré). — Un disque ouvert non vide con- 
tenu dans le disque unité est dit équilibré si et seulement si sa distance au cercle 
unité est égale à son rayon. Pour c G D{0, 1), on note De le disque équilibré de 
centre c. Sont rayon est r = ^-j^. 

Si de plus \c\ > -g, posons c' = (|c| — §) x ]f[ et soit Dc> le disque équilibré de 
centre c' . Alors |c| — | = '''^'^^ > donc le rayon r' de D^' vérifie r' = jr. De plus 
1 — |c'| = |(1 — |c|). On dit que De' est le fils de De- Si \c\ < ^ alors le fils de De 
est par définition Dq = D{Q,^). 



JACOBIENS, JACOBIENNES ET STABILITÉ NUMÉRIQUE 



23 




Figure 1. Fils d'un disque équilibré 

Soit f{z) une série entière d'une variable et d'ordre de grandeur {A, n) avec A>1 
et n > 1. Soit De C -D(0, 1) un disque équilibre où / est majorée en module par 
< e < 2~^°°. On note r le rayon de D^. Soit D^' le disque fils de et r' son 
rayon. Le disque de centre d et de rayon r/2, est contenu dans D^- Donc / y est 
majorée en module par e. 

La formule de Cauchy majore les dérivées de / en c/. 



|/W(c')| 



k\ 



fie' + 



< e- 



On veut majorer |/| sur Do'- On choisit un entier positif u et on majore 
séparément la partie principale d'ordre u en c' notée Pc'.u et le reste Rc',u — f—Pc',u- 
D'une part 



0<k<u-l 



k\ 



(5) 



< 



0<fc<«-l 



2^ T TF^Mô 



D'autre part, posant z = c'+y{l — \c'\), le reste Rc\u{c' +y{l — \c' \)) n'est autre que 
le reste d'ordre m en de la série recentrée y i-^ Fc'{y). Puisque z = c' + y{l — |c'|) 
appartient au disque équilibré De', le vecteur y parcourt le disque équilibré Dq. 
En d'autres termes, \y\ < |. On applique les lemmes 16 de recentrage et 17 de 
majoration du reste. 

La série recentrée Fc' est d'ordre de grandeur (Ac',n+ 1) avec Ac' majorée par 

Suivant les notations du lemme 17 on pose 



Br/ = 



2(2e)"+2(n+ l)!n!A 



y.n+2 



< Aexp(e)n^(l + |logr|)). 



24 



JEAN-MARC COUVEIGNES 



Pour z e r>c' on a 

Soit alors k > le réel tel que log k = '12 Lg'^'^ ^ P^^^ petit entier plus 

grand que = 
On suppose que 

llog/î] > IlogSc'l 

donc u > ^(1 ^"^'^j^^^"^ ^e'I) ^ suppose en outre que 

16(n+l)=^ 
(log 2)2 • 

Alors 



On déduit de (5) que 

5 

log |Pc',«(c' + 2)1 < ulog- +loge. 

On a u < 1 + ;Li2|4 donc ulog | < log | + < car | loge] > 1001og2. 

Donc log|Pe',„(c' + zJ| < 

Ainsi log|/| <log(2max(|Pc.,„|,|iîc',„|)) < ^^(j^lf^ loge + log 2 < 0.05 loge 
car loge < —100 log 2. 

Lemme 18 (Prolongement au disque fils). — Il existe une constante positive 
effective ci^ telle que l'énoncé suivant soit vrai : 

Soit f{z) est une série entière d'une variable et d'ordre de grandeur {A,n) avec 
A> 1 et n > 1. Soit D = C £'(0, 1) un disque équilibré de centre c et de rayon 
r où f est majorée en module par un Q < e < 1. Soit D^i le disque fils de De- On 
suppose que — loge est minoré par ci2(log^ + n^ \ logrj). Alors f est majorée en 
module sur le disque fils De' par e^ô . 

A. 3. Polygone de Newton d'une série entière. — Soit / une série entière 
d'une variable complexe et soit R< 00 son rayon de convergence, supposé non nul. 
Soit r un réel positif inférieur à R. On note D = D{0, r). 

On s'intéresse aux zéros de / dans D. Combien sont ils? Où sont ils? Comment 
sont ils affectés par une petite perturbation de / ? 

Pour répondre à ces questions, on cherche à enfermer les zéros de / dans une 
collection finie de petits disques disjoints tels que le nombre de zéros de / dans 
chaque disque ne soit pas modifié par une petite perturbation. 

On étudie d'abord la situation autour de zéro. On suppose que /(O) = 1 donc 
f{z) = 1 + X]fc>i /fc^*^- note d le degré de / en z, qui est en général infini. Le 
nuage de Newton associé à / est l'ensemble de points (/c, — log |/fe|) pour k > et 
fk ^ 0. Le polygone de Newton de / est la fonction M àe [0, ci] dans M définie comme 
le maximum des fonctions affines : [0, fi] — > R qui passent en dessous du nuage de 
Newton (c'est-à-dire (t>{k) < — log|/fe| pour tout k). La fonction A/" est bien définie 
car le rayon de convergence R est non nul. C'est une fonction convexe de [0, d]. Elle 
est affine sur tout intervalle ouvert délimité par deux entiers consécutifs. 
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En effet, soit A; > un entier inférieur à d. Pour tout e > il existe une fonction 
afRne qui passe sous le nuage de Newton et telle que Af{k) — e < (j){k) < N{k). 
De même il existe une fonction affine V qui passe sous le nuage de Newton et telle 
que M{k + 1) — e < 'tl){k + 1) < M{k + 1). On définit la fonction affine de la 
façon suivante. Si 0(fc + 1) < ip{k + 1) et tp^k) < (j>{k) alors est la fonction affine 
qui vaut (j){k) en k et ^p{k + 1) en k + 1. Si (/>(A; + 1) > tp{k + 1) alors Kc = (/>. Si 
(j){k + 1) < 'tp{k + 1) et ilj{k) > 4>{k) alors = ip. On vérifie que Ke passe sous le 
nuage de Newton. Quand e tend vers la famille des Ke converge simplement sur 
le segment [k, k + 1] vers la fonction affine ko qui vaut J\f{k) en k et J\f{k + 1) en 
k+1. Donc Ko minore J\f sur cet intervalle. Un argument de convexité montre qu'on 
a égalité. 

Ainsi est continue et affine par morceaux sur [0, d]. 

Les sommets du polygone de Newton sont les discontinuités de M' plus (0, 0) et 
éventuellement {d, J\f{d) ) . 

Soit k un entier entre et d. On pose l = N'{k). On appelle tangente en {k,l) 
au polygone de Newton, toute droite passant par P = (fc, l) et qui passe sous le 
nuage de Newton. On note a~ la dérivée à gauche, qui est la pente de la tangente à 
gauche. Le vecteur (— a~, 1) est orthogonal à cette droite et tourné vers l'intérieur 
de Af. De même a+ est la dérivée à droite. On suppose que a"*" > a~ donc P est 
un sommet. Soit a dans ]a~, «"*"[. La position relative du nuage de Newton et de 
la tangente en P de pente a nous renseigne sur l'ordre de grandeur de f{z) pour 
un z tel que log \z\ = a. En effet pour tout tel z et pour tout entier positif m on a 
-logl/m^™! = -logl/^l -ma = (-a, 1) • (m, -log 1/^1) > (-a, 1) • (fc, - log 
De sorte que /feZ*' est le terme dominant sur le cercle \z\ = exp(a). Il reste à voir 
jusqu'à quel point. On se doute que si le sommet P est assez anguleux, les autres 
termes peuvent être négligeables. 

Soit donc m ^ fc un entier positif ou nul. Le point (m, — log |/m |) est au dessus 
du polygone de Newton. Si m > fc il est donc au dessus de la droite passant par 
P et de pente a"*". Donc \ fm\ < |/fe| exp(— (m — k)a~^). Donc, pour log|2;| = a, le 
terme fmZ™ est majoré en module par |/fc||2;|'^ fois exp((m— fc)(a — a+)). La somme 
X]m>fe l/felkl*^ donc majorée par l/feU^j'^ fois en posant x = exp(a — q:+). Si 
m < k le point (m, — log \fm\) est au dessus de la droite passant par P et de pente 
a~. Donc \ fm\ < \fk\ exp(— (m — k)a~). Donc, pour log \z\ = a, le terme fmZ"^ est 
majoré en module par fois exp((m — k){a — a~)). La somme X^^^fc |/fc||2:|*' 

est donc majorée par l/feUzl*^ fois en posant y = cxp(q;~ — a). 

Pour toute fonction h holomorphe sur un voisinage du disque fermé de centre et 
de rayon exp(a), et non nulle sur le bord de ce disque, le nombre de zéros de h dans 
l'intérieur de ce disque est l'indice en du lacet image du bord {z, \z\ = exp(a)} 
par h. C'est un invariant homotopique discret. Donc dans toute famille continue de 
fonctions holomorphes sur un voisinage de ce disque et jamais nulles sur son bord, 
le nombre de zéros dans l'intérieur du disque est constant. 

On en déduit que si 73- + -r^ < 1 le nombre de zéros de / dans le disque ouvert 
L'(0,cxp(a)) = {z, \z\ < exp(Q;)} est le même que celui de z'' soit k zéros. Cette 
condition est satisfaite si a; et y sont inférieurs à ^. 

On appelle pente du polygone de Newton, une valeur de la dérivée de J\f en un 
point ou elle est dérivable. 
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Lemme 19 (Polygone de Newton). — Soit / = 1 + X]fc>i /fc-^*^ ^''^^ série 
entière de rayon de convergence R> 0. Soit ^ G D{0, R) un zéro de f. Il existe une 
pente a du polygone de Newton telle que \ log |Ç| — a\ < logS. 

On note V3 l'intervale ] — oo,logiî[ privé des intervales [a — log 3, cr + log 3] où 
a parcourt l'ensemble des pentes du polygone de Newton. 

Si a est un réel de P3 il existe un unique sommet P = {k,J\f{k)) admettant 
une tangente de pente a. La fonction f a exactement k zéros dans le disque ouvert 
£>(0,exp(a)). 

On note V4 l'intervale ] — 00, log /?[ privé des intervales ]a — log4, a + log4[ où 
a parcourt l'ensemble des pentes du polygone de Newton. 

Si a est un réel de V4 il existe un unique sommet P = {k,Af{k)) admettant une 
tangente de pente a. Pour log \z\ = a on a \ f{z)\ > = SSltîïi. 

Comme on pouvait s'y attendre, ce lemme est moins précis que son pendant non- 
archimédien. On ne peut pas l'utiliser directement si les pentes sont trop proches 

les unes des autres. Dans ce cas, on pourra former (par exemple) la série g(z) = 
î{\fz)f{—\fz) dont les zéros sont les carrés des zéros de /. Le passage de f k g 
clarifie la situation dans le voisinage du cercle unité. On peut réitérer l'opération si 
nécessaire. 

A. 4. Le plus petit zéro d'une série entière. — Soit F = Fq + J2k>i ^kz'^ 
une série non constante de rayon au moins 1 telle que Fq ^0. La série normalisée 
/ = F/Fo admet au moins une pente. Soit tri la plus petite des pentes. On suppose 
d'abord que cri est négatif. Soit alors logr la borne inférieure de ](Ti,0[n7'3. Si 
ce dernier ensemble est vide on pose r = 1. Si r < 1, alors / admet un zéro de 
module < r. On veut montrer que si Fq est petit alors r est petit ou bien F est 
uniformément petite. On suppose que |Fo| < 1. Le segment ](Ti, logr[ est couvert par 
des intervales fermés de rayon log 3 centres en les pentes du polygone de Newton. On 
note CTi < (72 < • • • les pentes successives. On a a2 < CTi + 2 log 3, < ai +4 log 3, 
. . ., cT/s < C7i + 2{k — 1) log 3, tant que ai + {2k — 3) log 3 < logr. On pose donc 

_ logr - 0-1 +log3 
' 21og3 ' 

et pour tout 1 < k < £ on a, ak < ai + 2{k — 1) log 3 et donc J\f{k) < kai + k{k — 
1) log 3. Cela prouve en particulier que le degré d de F est au moins égal à L 
On pose k = l ei on obtient 



(6) A/'(^) < ^cri l)log3. 

Le principe de la démonstration est le suivant : on suppose Fq petit. 

Si ai est grand alors la série F est petite car ses premiers coefficients sont petits. 

Si ai est petit et £ petit alors r est petit : le polygone de Newton est anguleux 
près de l'origine et il y a une petite racine. 

Si ai est petit et £ grand alors la pente du polygone de Newton varie peu au 
début, et la série F a de grand coefficients. 

Pour formaliser ce raisonnement, nous supposons maintenant que F est d'ordre 
de grandeur (A, n) avec A > 1 et n > 1. Donc — log \fk\ = — log + log |i^ol est 
minorée par — log A — nlog(A; + 1) + log |Fo| qui est une fonction convexe de k et 
qui minore donc le polygone de Newton. Pour k = t on obtient 
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(7) > - log A - n log(f + 1) + log |Fo I . 

■-<Ti+3 
21og3 



Comme l < ^"sr-a-i+aioga a cti < -2£log3 + logr + 31og3. En reportant dans 



l'inéquation (6) on a Af{£) < — ^^logS + £(logr + 21og3). L'inéquation (7) donne 
alors£2log3-£(logr + 21og3) < log A + nlog(^+ 1) -log |Fo| < log A + n£-log |Fo|. 
Donc £ satisfait l'inégalité quadratique 

(8) f^log3-^(logr + 21og3 + n) -logA + log|Fo| < 0. 

Si £ > logr + 2Iog3 + n alors on déduit de l'inéquation (8) que £^(log3 — 1) < 
log A - log |Fo I • Au total 



(9) £ < ma.(logr + 21og3 + n, f-^^^^^^^)- 

Comme £ > '°^7iogt'°^^ on a en > -2^ log 3 + logr + log 3. On déduit de 
l'inéquation (9) 



<Ti >min((l-2 log 3) log r-2n log 3-(2 log 3) = +log 3,-2 log 3^ +log r+log 3). 

On rappelle que r < 1, et on suppose que 

-log|Fo| > 0(logA + n^). 

On en déduit alors que 



cTi > -lOV- log |Fo| + logr 



Si — logr > \/— log |Fo| on s'estime heureux puisqu'on a montré que F admet 
un zéro très petit. Sinon on a cri > — ll-\/— log \Fq\. On observe que cette dernière 
inégalité est vraie aussi si ci est positif ou nul. 

Donc pour tout entier k > on a fk = < cxp(— fccri) < cxp(llA:-\/— log |Fo|). 
Si z € D{0,^) est un complexe de module inférieur à ^ alors iF^Hzl*^ < 
|Fo| exp(/;;(ll-y^^IôgjF^ — log 2)) et pour tout entier positif u, la partie principale 
Pu{z) = Y,o<k<u majorée en module par |Fo| exp(u(ll-^- log |Fo|-log2)). 

Si on choisit u = alors log|P„(z)| < log|Fo| +u(-log2 + 

lV-log|Fo|) < iHim. 

Pendant ce temps là, on peut majorer le reste Ru{z) à l'aide du lemme 17. On 
demande que u > (\'^^2)^ '^^ acquis si 

-\og\Fo\>On\ 

Suivant les notations du lemme 17 on pose B = 2"+^n!v4. On peut majorer Ru{z) 
en module par k pourvu que — log k < — log B+u^^^. On note que u = 2if^^ ~J 

est plus grand que 23^^" " \/— log \Fq\ > O. 
Si v'-logl-Fbl > ©(log A + V?) alors 



logB< i^ v/-log|fb| 



donc - logB + > V-logl^^ol X 
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On pose donc — logK = 0.007-\/— log |Fo|. 

Comme la partie principale Pu{z) est majorée par -y/lFol qui est plus petit que k 
on alog|F(z)| < log 2 - 0.007a/- log |Fo| < -0.006a/- log |^b| si ^~\og\Fo\ > O. 

Lemme 20 (Plus petit zéro). — Il existe une constante positive effective cis 
telle que l'énoncé suivant soit vrai : 

Soit F{z) = Fq + Fiz + ■ ■■ une série entière d'ordre de grandeur {A,n) avec 

A > 1 et n > 1. On suppose que \Fa\ < 1 et log |Fo| > Ci3(n^ + log A). Alors 
ou bien f a un zéro ^ tel que log|Ç| < —a/— log |Fo|, ou bien f{z) est majorée en 

module pour z € £'(0, 5) par k tel que logn = ^ 200^ 

A. 5. Stabilité des zéros d'une série entière. — On se donne maintenant une 
série entière / d'ordre de grandeur (A, n) avec A > 1 et n > 1 et un complexe c tel 
que |c| < 1. On veut montrer que si /(c) est petite alors c est proche d'un zéro de 
/ ou bien / est petite sur le disque i?(0, |). 

On note F^iy) = /(c + y(\ — |c|)) = Fq + F\y + • • • la série recentrée en c. Elle 
est d'ordre de grandeur (Ac,n + 1) avec A^ = ^n!exp(n + 1)2"+^(1 — |c|)~"~^. 
Appliquons le lemme 20 à la série recentrée Fc- On suppose donc que Fq = f{c) 
vérifie |/(c)| < 1 et 



V-log|/(c)| > ci3(n2 + log 



Il vient que / a un zéro ^ tel que log|c — Ç| < —\/~ log |/(c)| ou bien / est 

majorée en module par /t sur le disque équilibré où log k = — ^ ^200 ^ ~ • -D^'^s 
ce dernier cas, on peut appliquer le lemme 18. Soit w = [— ^^^^Ç^l et soit v tel 
que log V = et supposons que 

-logz/>ci2(logA + n2|log( i^-^ ) I). 



2 

On applique w fois le lemme 18 et on montre que / est majorée en module par 
V sur le disque £)(0, |). 

log 20 

On note que 20"' < 20(1 - |c|) < 20(1 - |c|)-". 



Lemme 21 (Stabilité d'un zéro). — // existe une constante effective positive 
Ci4 telle que l'énoncé suivant soit vrai : 

Soit f une série entière d'ordre de grandeur {A,n) avec A> 1 et n > 1. Soit c 
dans £>(0, 1). On suppose que e = |/(c)| < 1. Soit v tel que \o%v = 

On suppose que — logv > ci4(log A + n^(l + | log(l — |c|)|)). 

Alors f admet un zéro £, tel que log |c— ^| < —\/— log |/(c)| ou bien f est majorée 
en module par v sur le disque D(0, |). 

Puisque les zéros de / bougent peu sous l'effet d'une petite perturbation, ils ne 
doivent pas beaucoup s'éloigner des zéros de la partie principale. Et c'est un moyen 
commode de les localiser. 

Soit donc / une série entière d'ordre de grandeur (^, n) avec ^ > 1 et n > 1. 
Soit r un réel entre 1/2 et 1. Soit e un réel entre et 1 et soit v im réel positif 
tel que logj^ = — ei'^~'°^'^ avec C15 constante plus grande que 1 et C14. On 
suppose que |/| n'est pas majorée par v ^ e sur -D(0, i). On suppose que — log:^ > 
ci5(log^ + n2(l + I log(l - r)|)). 
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On cherche un entier positif u tel que le reste Ru{z) soit majore en module par 
e sur L>(0, r). Selon le lemme 17 il faut que u > avec B = j^^^p^- Si 

Ci5 est assez grand alors 

I logB| < I loge] 

donc il suffit que u > Soit donc u le plus petit entier satisfaisant cette 

condition. On vérifie que u > '^11°^^! > (log" )^ '^is assez grande. Donc Ru{z) 
est majoré par e en module sur D(0, r). On pose p = exp(— loge). 

La partie principale P„(z) est un polynôme de degré m — 1 qui a donc u—1 zéros 
dans C. Donc l'intervalle [r — 4up, r] contient^*) au moins un réel positif R tel que 
\R, — ICII > 2/9 pour tout zéro ^ de Pu- Soit V le fermé de C obtenu en retirant au 
disque fermé D{0,R) tous les disques ouverts D{£^,2p) où les ^ sont les zéros de 
Pu{z). D'après le lemme 21 le polynôme Pu{z) est strictement minoré en module 
par e sur le domaine formé V. Donc f{z) et Pu{z) n'ont pas de zéro dans T). Elles 
ont le même nombre de zéros dans -D(0, R). Elles ont le même nombre de zéros dans 
chaque D{^,2p). Donc les zéros de Pu{z) dans D{0,R) approchent ceux de f{z) à 
distance 4p. On obtient le 

Lemme 22 (Stabilité globale). — Il existe une constante effective positive cie 
telle que l'énoncé suivant soit vrai : 

Soit f une série entière d'ordre de grandeur {A,n) avec A> l et n> l. Soit r 
et p deux réels tels que ^ < r < 1 et < p < 1. Soit u le plus petit entier plus grand 
que "^ij"!^! . On suppose que —(1 — ry^logp > ci6(logA + n-^| log(l — r)|)). Alors 
f satisfait l'une au moins des deux propriétés suivantes : 

1. Sur le disque £)(0, 5), le logarithme log|/| du module de f est majoré par 

C16 

2. Il existe un réel positif R tel que r — Aup < R < r et tel que dans le disque 
D(0, R) les zéros de f{z) sont approchés à distance 4p par ceux de la partie 
principale Pu{z) de degré u—1. En particulier, il y a au plus u tels zéros. 

Théorème 3 (Zéros d'une série). — Soit f = J2r>ofr^^ série entière 
d'ordre de grandeur {A, n) avec A>letn>l. Soit R un réel compris strictement 
entre et 1. Soit fi la partie positive de — logmax|^|<i |/(2;)|. Le nombre de zéros 

de rn,odule < R est polynomial en n, log^, (1 — R)^^ et /;,. 

// existe un algorithme qui pour f et R comme ci-dessus^^^ et pour k entier positif , 
retourne 

- un rationnel R' tel que \R — R'\ < exp(— fc), 

- le nombre de zéros de f dans D{0,R'), 

- une approximation de ces zéros à exp(— fc) près, 

en temps déterministe polynomial en n, logA, {1 — R)~^ , p et la précision absolue 
k requise. 

Cela découle du lemme 22. Il suffit de rappeler l'existence de tels algorithmes 
pour la recherche des racines d'un polynôme. □ 



(*)Si C15 est assez grand alors 4iip est (beaucoup) plus petit que r, donc l'intervalle en question 

est constitue de réels positifs. 

(^'La série / est donnée sous la forme d'un oracle qui calcule les coefficients fr en temps polynomial 
en r et en la précision absolue requise. 
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